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The combination of high frequency vibrations and metalhc transport in graphene makes it a 
unique material for nano-electromechanical devices. In this letter, we show that graphene-based 
nano-electromechanical devices are extremely well suited for charge pumping, due to the sensitivity 
of its transport coefficients to perturbations in electrostatic potential and mechanical deformations, 
with the potential for novel small scale devices with useful applications. 
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Device miniaturization has led to small size mechanical systems, NanoElectroMechanical devices (NEMs) with a 
wide range of uses in fundamental and applied research [IH5]. In particular, electron pumps and turnstiles have been 

— I extensively studied[4 6J, including NEMs based devices [7inT|. Graphene NEMs fT^HT^ have an enhanced tunability 

—i with respect to devices based on carbon nanotubes, while keeping advantageous features such as high vibration 

2 frequencies and metallicity. Suspended graphene samples have a very high electron mobility[T3], and a large and 
well characterized electronic coupling to the strains induced by long wavelength vibrations [16j. Long wavelength 

C^ strains in a ballistic graphene sheet modify the electronic transport coefficients through the sheet [T7]. A flexural 

^ deformation leads to uniaxial strains within the suspended area, inducing a strain mismatch at the boundary between 

tn the suspended and non suspended regions, modulating the transport coefficients. Deformations of amplitudes of 

.4_i a few nanometers in samples of microns in size and the tuning of its electrostatic doping can be simultaneously 

2 achieved by adjusting the electrostatic force between the graphene layer and the metallic gate below it 17;. The 

fij periodic modulation in time of these internal parameters, i.e. electrostatic doping and strains, make possible to 

'■ achieve adiabatic charge pumping |18H20j . if the appropriate symmetries are broken. We argue below that these re- 

pH quirements can be met in realistic experimental setup, leading to charge pumping of the order of few electrons per cycle. 

O 

CJ We analyze the feasibility of a pumping device using the geometry sketched in Fig. [I^. The length of the sheet 

'~~' is L, and the applied voltage is V{t) = Vdc + VacCOs{u)t). We describe the deformation in terms of a single degree 

I of freedom, the maximum vertical displacement, a{t). Its dynamics is determined by the sum of the time dependent 

^ electrostatic force between the sheet and the gate, J-e, the restoring elastic force, J-s, and a dissipative term introduced 

00 phenomenologically, .FnPT]: 

^ 5'« X- ^ ^ ^ X- 

O ^9^2" = J-S+~f-D+J-E 

Cn ^ C^VdcVac , ,. 

^-5 J-E = COS(wt) 



^s - - y ^4 "^ (a^ + ia^ho + 3a/tg) + -^ (A + 2^)a 

d where p is the mass density, A and /i are Lame elastic constants, Ct is the total effective capacitance due to the 
back-gate oxide and air dielectric, AL and h^ describe the amount of slack and vertical displacement of the sheet in 
the absence of the periodic driving potential. The phenomenological parameter r^ describes damping, and the quality 
factor is Q = (woTd)/2, where wq is the resonant frequency. Gurrently, experimentally obtained wq for graphene is in 
the range of 100 MHz [T31II11IH]. Fig. [lb reproduces a typical experimental ujq as function of Vdc with our model. In 
the linear response regime, wq « ht^jlrp^l'^ , whereas ft-o can be tuned through Vdc and is proportional to {r?L'^Y^'^ ■ 
Gontinual device downscaling and improvements in graphene fabrication processes will allow for GHz operation. 
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already realized in nanotube systems [25], 

We look for the frequency and phase response to the dynamical system described by Eq. [T] The equations define 
a non-linear resonator, which we solve approximately ;21 using techniques derived for the Duffing modei p^ I25j . We 
show in Fig. [Tp the dependence of the maximum amplitude, a{uj), for different driving force Vac- When the driving 
force exceeds a given threshold, the oscillator shows bistability and hysteresis jT4]. Our results are in reasonable 
agreement with experimental data[14 shown in inset. Time varying deformation of graphene modifies its electronic 
spectrum through the modulation of electrostatic doping and in-plane strain modeled with, 

1 

Uxx{t) = U2;a;{l + <^'"xxSin(wt + (/))} — (2) 

where Sdg is the Dirac point energy in graphene with Fermi energy taken as zero, and e^ = hv ^{TiCTVdc/ c)^ , 
Std = Vac/Vdc, Uxx — SHq/SL^ and Su^x — a/ho. The internal parameters, £d and Uxx, constitute the two parameters 
for adiabatic quantum pumping in graphene NEMs, and are governed by the amplitude and phase response of the 
resonator system. Fig. [ip-e shows the dependence of amplitude a{uJo) and the phase response 0(a;o) on Vac and 
the quality factor Q. Improvements in quality factor, where values as high as Q = 10^ at T = 90 mK have been 
reported|26j. will lead to stronger non-linearity and sensitivity. 

Cyclic variation of the two internal parameters given by Eq. [2] constitute a scheme for quantum pumping. The 
scattering wave ipj^x) in the various regions: left contact, graphene and right contact, denoted by the subscript 
j — £, g,r respectively, can be written as follows: 



i^j{x) = < 



1 pik^ex I -7? ( ,1 p-ik.j;tx 

ai I I e*'^="9^ + cko I +1 e~*'^="9^ (3) 






Here, r]j are the pseudospin phases defined as, rjj = Hvf l^^g"^ where £dj is the Dirac energy in each region. TZy, 
Tv, cti and ag are the wave amplitude coefficients, to be determined by imposing wave continuity at the interfaces. 
The in-plane strain Uxx leads to an effective gauge potential[17|. Ay = zLng ^^^^*" where /3 = —-g^(§) ~ 2, i^ ~ 3eV 

is the nearest neighbor hopping term, few 1.4 A is the bond length, n^ is a dimensionless geometrical factor which is 
found numerically to be ~ 0.5, and the two signs correspond to the two inequivalent Dirac points in the Brillouin zone 
i.e. K and K' . It modifies the transverse wave- vector through hkyg = hky — eAy. Time varying transport coefficients 
TZy{t) and %(t) are determined adiabatically from Eq. Is] The pumping current for each valley is |20l[28] . 

^^ = ^SEy ^v d6^i7„(fc„i) (4) 

where v denotes the valleys (i.e. K,K'), /o(e) is the Fermi-Dirac distribution and the pumping coefficient is defined 
as, fly — '^T^ + '^^TZl Evanescent contributions, albeit small, are also included in the model. 

In order for the pumping current to be non-zero, spatial inversion symmetry needs to be broken. Typical charge 
pumping scheme employs two electrostatic gates to achieve this[29]. In NEM-based quantum pump, a number of 
perturbations will achieve that. In the following, we assume that the left and right contacts are not equivalent, which 
is modeled by different densities of states. In reality, this can be implemented by using different materials for the two 
contacts(3D]. We assume ballistic transport, which implies that the mean free path, £, is larger than the dimensions 
of the device, i > L. This limit can be achieved in clean suspended samples[l6]. Diffusive scattering will suppress 
the effect of the gauge field 1 17). so that the modulation of the scattering matrix will be reduced, but, for sufficiently 
low amounts of disorder, a finite pumping current will exist. 

Using the model presented above, we consider a prototypical device of L = 50 nm, AL = nm and W = 1 /im. 
Symmetry of the problem requires that the Hamiltonian HKiky) = T-LK'i—ky) {y is aligned along the zigzag 
direction), which also implies TK,k = ^K',-k ■ In other words, the pumping current ly = ^j, ly^k from valley 



V = K, K' must be equal and flows in the same direction. Hence, in subsequent analysis, we shall consider only 
one of the valleys i.e. K. First, we illustrate some of the basic features of electron pumping in graphene NEMs. 
Fig. [2k-b plots the transmission Txi^y) and pumping coefficient ilpciky) over a pumping cycle for = (top panels) 
and 9 = 7r/2 (bottom panels). In these calculations, we assumed an asymmetric contact doping of £m = —OAeV 
and Sdr — ^O.SeV. The contact with a lower doping will stipulate the maximum allowable transverse momentum 
wave- vector (kmax) that could accommodate propagating states through the device. As the graphene resonator 
undergoes strain modulation, it induces a translation in its transverse momentum hkyg = hky — eAy. States where 
kyg > kmax would be evanescent in the contacts and their transport coefficients will be zero i.e. white regions in 
Fig. [2k-b. In general, larger ky states leads to stronger interference effects as seen in Fig. [2k. Since pumping current 
is proportional to the accumulated complex phase per cycle, Qk is most significant at larger ky. When the two 
parameters are in phase, ^k for a given ky state is exactly antisymmetric within each time cycle, i.e. the ^ — >■ ^ is 
anti-symmetric with —§—>■§ portion of the cycle, hence Tk = 0. This symmetry is broken when 0^0, and a finite 
pump current then ensues. 

Fig. [2t-d plots the time averaged conductance (G) and the pumped charge per cycle Qc for varying transverse 
momentum, ky, and doping, £dg- Here, we observe a larger flu at negative ky and vice versa for K' valley i.e. a 
valley Hall effect. Based on the condition iK.k — 1-K' -k stated earlier, it is apparent that a valley Hall effect 
will be present, since T^.k 7^ ^K.-k in general. The valley Hall effect will induce a spatially dependent valley 
polarized current, whose effect is maximal near the two edges. Calculations as shown in Fig. [2tl estimate the valley 
polarization, i.e. {Ik — Ik')/Ik, to be as large as 90%. Fig. [2^-g show that the pumped charge Qc is finear with 
respect to the amplitudes of the pumping parameters and the device length. The latter is a result of increasing 
interferences frequency with L. Q^ also increases with contacts doping asymmetry, except that the effect maximizes 
when density-of-states in one of the contacts becomes the bottleneck to conduction. Reasonable driving voltages lead 
to measurable currents for devices with similar features to experimentally studied NEMs. These systems provide a 
robust setup where quantum pumping can be observed. 

We brieffy discuss issues related to experimental realization. In conventional quantum pumping scheme, displace- 
ment current induces by stray capacitances can interfere with the quantum pumping dc current [S^ I31j . as the two 
gates can work in unison to result in a rectification of the displacement currents 32 . Since our proposal utilizes only 
a single back gate, there will be no rectification of the ac displacement currents at least to first order in frequency. 
The calculated values of the current in our device are such that situations where the charge pumping per cycle is 
close to one or a few electrons are feasible. Coulomb blockade effects will favor the transference of an integer number 
of electrons per cycle, so that the ratio between current and frequency will be quantized. Such behavior will manifest 
itself as steps in the dependence of this ratio on driving voltage. The charging energy of a device of length L is 
Ec « e^/L, so that Ec ~ lOK for L ^ l/^m, and Coulomb blockade effects can be expected to be relevant at lower 
temperatures. The observation of quantized steps in I/u) will allow for the realization of a graphene based current 
standard [33 , making graphene an unique material from whom current and resistance [34J standards can be fabricated. 
Note also that the carrier density in very clean suspended graphene samples can be adjusted with great accuracy, 
making the physics at the Dirac point accesible|T2]. At these concentrations, electronic transport in ballistic systems 
is determined by evanescent waves [35l I36j. and pumping through these modes can also be expected [3 7]. In principle, 
we also envision alternative schemes via optical means [38 , where the laser could induces a non-equilibrium electronic 
temperature which through coupling with the flexural phonons will lead to strains and vibrations. 

In summary, we show that a graphene NEM near resonance can function as an adiabatic quantum pump under 
realistic experimental condition, due to the unique electronic coupling to the strains induced by long wavelength 
vibrations. Experimental realization of this effect would open up new opportunities in fundamental and applied 
research with graphene NEMs[I]-E]- 
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FIG. 1: (a) Schematic of a typical graphene nanoelectromechanical resonator actuated electrostatically with a back gate. Gating 
capacitance is given by the total effective capacitance due to the back-gate oxide and air dielectric i.e. Ct = [^(7 ('i + '^o) + e5i02*S«02]~^i 
where we assumed tsi02 = 200 nm and d = 100 nm in this work, (b) Resonant frequency /o as function of bias voltage Vdc, computed using 
our model i.e. iJo = y^kojp, where kg = daJ^s{a = 0) is the linearized spring constant term, p and AL are used as fitting parameter to the 
experimental data of 2 devices (in symbols) reproduced from 1 14). (c) Amplitude response, a{uj), of device 1 for different driving forces Vac, 
obtained by solving the non-linear resonator model of Eq. IT] using techniques employed for the Duffing model, assuming a quality factor 
Q = 125, the value corresponding to the experimental situation|14|. The oscillator shows features of bistability and hysteresis similar to 
that of experiments 1 14| (see inset and Ref. 1141 for measurement details), (d — e) Amplitude and phase response at resonance (of device 
1) as function of driving force Vac for 2 different quality factor Q = 25 and 125. 




(a) Traiisnussioii T (b) Piiinpuig coeflicieiit D 

yT imx \ 




0.0 0.2 0.4 0.6 -0.4 -0.3 -0.2 50 100 150 
Suxx Ej. (eV) L (mn) 

FIG. 2: We consider graphene NEMs based electron pumping device, through cyclic variations of £.^(1) and XAx-xify as described in 
Eq. [2] Unless stated otherwise, we consider graphene dimension of L = 50 nm, AL = Onm and W = 1 ^m, with equilibrium parameters 
Uxx = 0.02 and e^ = — 0.2eV. Contact asymmetry is introduced through S^ = —OAeV and £,1^ = —O.SeV. (a) Transmission, Tfc{ky), 
as function of time over one pumping cycle, for cases where the two parametric variations are in-phase (i.e. </< = 0) and out-of-phase (i.e. 
(p = it/2). In these calculations, we assumed Suxx = 0.8 and SEg^g = 0.2. Dashed lines indicate the minimum and maximum transverse 
momentum ky (black) and ky — iAy (red), (b) Similar to (a), except for pumping coefficient Qidky)- Note that pumping current for 
the 4> = case is zero, (c — d) Time averaged conductance (G) and charges per cycle Qc as function of graphene's doping ^'^jg and 



transverse momentum k„ 



Dashed lines indicate ±hvf{ky — i^y)- 



In these calculations, we assumed 5ux 



0.2 and SS^g = 0.2. Note 



that calculations for (a-d) arc performed for only one of the valley i.e. K. (e — g) studies Qc as function of various parameters: pumping 
amplitude Suxx, contact doping asymmetry and device length L. In these calculations, we assumed Suxx = 0.2 and SS^ig = 0.2, unless 
stated otherwise. 
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We consider a suspended graphene resonator as described in [1, 2]. The dynamical equation for the out-of-plane 
deformation from equiUbrium a{t) is given by (see the respective Appendixes on the form of the forces), 

p^ = ^s + :fd + :fe (1) 

where Tg, J-jj, Te are the time dependent restoring elastic, damping and electrostatic forces (or pressures) respectively, 
in units of to~^. p is the mass density. The total restoring elastic force, ignoring second order terms due to bending 
forces K and 0{h^), reduces to, 

Ts = -Y^T^ {a' + 3a^ho + 3ahl) + ^-^{X + 2p)a (2) 

obtained by minimizing elastic energy. Damping is treated phenomenologically via, 

Td Ot 

Td can be easily obtained from quality factor Q measured in experiments. Q is defined to be Q = wo/Aw = ujQTd/2, 
where Aw is the so called bandwidth of the resonance peak. Reported Q is around 100 [1, 2] at room temperature, 
depending on many factors. For example, temperature dependence of Q ex T^^-^^ was found, for T < 1{)QK [2]. And 
a record Q k, 100,000 at 9QmK was reported [3]. Mass density of graphene assumed to be p « 7.4 x IQ^^kgmT^ . 
Lastly, the electrostatic force is modelled as, 

Te = —^ cosiiot) (4) 

eo 

neglecting 0{V^^) and other non-linear terms. Ct ~ [f-o^{d + ho) + eg/o2^Si02]^^ is the total effective capacitance 
due to the back-gate oxide and air dielectric, d being the perpendicular distance of the unstrained graphene from the 
substrate. 

We are interested in the steady state solution to the dynamical equation i.e. a(a;) — \a\exp{i(j)) . We seek an 
approximate solution through an iterative technique used in the Duffing model [4, 5]. For convenience, we rewrite the 
dynamical equation as, 

pa = —k^a — fcifl — k2a a + fcos{ujt) (5) 

Td 

where 

6A{X + 2p)hl 8AL,. , , 6A{X + 2p)ho , 64(A + 2^^) , C^VdcVac ,^. 

ko = ^3 ^^(A + 2^) , k, = , k, = 3^4 ' / - -^^ (6) 

In the spring constant term ko, we consider only the case for AL < 0. k2 contributes to the Duffing force, and renders 
the spring more stiff (soft) if positive (negative). In the former, the effect will be a shift of resonance with increasing 
driving force. And at larger driving force would lead to bistability and hysteresis [2]. 



The frequency response a{uj) around the resonant frequency ujq — \/ko/p has the fohowing approximate solution, 



iu;l-u')\a\ + ^-^\a\ 



Ap 



iancj) : 



w|a| 

Td 

a;|a| 



Td 



H 
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Note that a^ terms affects the higher harmonics 2ujq. These converge to the Lorentz model solutions if we set k2 — 0. 
In the linear regime, i.e. ^2 = 0, the response at a; = cjq goes as \a\ = 2Qf / puj^. When k^ ^ 0, \a\ follows. 
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From this, we can define a threshold driving force /^/j where \a{u) ~ wq)] starts to deviate from linearity i.e. \a\ oc /, 
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As evident, p has no effect on ftt, and a larger Q yields a smaller fth i.e. increased sensitivity to non-linearity. 
However, a larger Q is desireable to acheive a larger resonance \a{uj = wo)|. 



APPENDIX A: ELECTROSTATIC FORCES 



The electrostatic force is modelled as. 
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neglecting Oiy"^^. Ct — [eo^('^+ ^o) + ^si02^Si02\ ^ is the total effective capacitance (per unit area) due to the 
back-gate oxide and air dielectric, d is the perpendicular distance of graphene from the substrate when unstrained. 



APPENDIX B: ELASTIC FORCES 



The position of a 2D membrane can be described by the in-plane and out-of-plane deformation field given by 
\i{x^y) — [ux{x,y),Uy{x,'y)\ and h{x,y). In the linear approximation, the strain tensor is given by. 



"a,3 = -^{daUfi + dpUa + dadph) 
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The elastic free energy is given by [6, 7 
£ ~ I dxdy 

dx 



-K(V^hf + -A(u,, + uyyf + A*(wL + «L + 20 - ^Eh 

1 ,_9,.9 1 



K{V^hf + -{\ + 2p)ul,- T%'h 



(B2) 



where k ~ leV is the bending rigidity, p ~ 9eVA^^ and A « 2eVA~''^ are the Lame constants of graphene. 
Experiments [8] measures an elastic constant for graphite cn — 106 x IQ^^NrnT^ « ITPa. For graphene, we have 
the relation cudin = A -I- 2p. Using an interlayer separation distance of dm = 0.335nm [9], we cudin « 355Nm~^. 
This is in good agreement with the values of Lame constants we assumed. Recent measurement [10] of the Young's 
modulus of graphene yields E2d = 3A2Nm~^ . Note that by definition, E2d = cudm- 

T^^ is the pressure induced by the bottom electrostatic gate. Assuming homogcnuity along transverse direction, 
we arrive to a one-dimensional problem. Minimizing £ in Eq. B2 yields us the following Euler-Lagrange equations, 

"3 
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where the latter imphes the longitudinal strain is constant i.e. 

Uxx = constant (B4) 

Finally, the differential equation governing h{x) is, 

Kd^h - (A + 2^i)uxxdlh ^ J-|,« (B5) 

The differential equation for h then reduces to, 

-{\ + 2ij)uxxdlh^F%^ (B6) 

whose explicit solution with the boundary condition h{±L/2) ~ is, 

8(A + 2/x)Ua;a: L^ 8(X + 2^)Uxx 

where /iq is the maximum deflection i.e. at x = 0. Including the second order terms, the expression can be rather 
complicated [11]. With the profile h{x), simple geometry allows us to relate Uxx with hg as per Eq. C9. Then the 
elastic force equation reduces from Eq. B6 to, 

•^l' = -y(A + 2A.)§ + ^(A + 2^)/.o=-Ff (B8) 

where we included the possibility of an initial tension, AL < 0. Elastic forces due to deformation a away from the 
equilibrium can then be described by, 

:Fs = -f ^^ {a' + 3a^h, + 3ahl) + ^(A + 2M)a (B9) 

APPENDIX C: DEFORMATION AND GATING 

Next we shall determine Uxx- Following [12], it is defined as, 

_ L' - {L + AL) V -L AL 
""" " lTAL ~ ~L L~ 

where L' is the length of the strained graphene, 

f-i/2 i'L/2 ^ _L3Cjre9^2 
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L' = 2 dx^lTWW^L+ dx{dxhf^ )^> ^ +L (C2) 

L is trench length and L + AL is length in absense of strain. Hence, solving Uxx then reduces to finding the root of. 



U^^ H ^u^^ — ::rTT-, :^^^ = (C3) 
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No initial tension or slack: If AL = 0, wc will obtain. 
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Electrostatically, one can approximate J^p as 



n' « -^ (C5) 
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where n is the carrier density in graphene and e^ff is the effective dielectric due to air gap and back gate oxide. 
We ignores the curvature of graphene. Assuming that AL — 0, we get an expression for maximum deflection ho at 
equilibrium, 
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With initial tension or slack: For general case of AL ^ 0, To obtain /iq rigorously requires solving the following 
the electrostatic and elasticity equations self-consistently. 
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When F'^ = 0, the latter requires ho = for AL < (tension), which is expected. Another remark. A series 
capacitance Cox to Ca is essential as it provides stability to the electrically actuated mechanical system. A relation 
between u^x and ho can be obtained from Eq. B7 and C3, yielding. 
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